Signatures of the BCS to Bose crossover in atom shot noise correlations 
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We discuss the superfluid wave function of an atomic Fermi gas, which can be measured in an 
atom shot noise correlation experiment. We note that it is qualitatively different in the BCS and 
Bose limits and on the basis of this we propose quantitative criteria to identify the crossover between 
the two regimes. 



The achievement of quantum degeneracy in Fermi 
ases 0, loading of ultracold atoms on optical lattices 
and control of atom-atom interactions using Feshbach 
resonances suggest using ultracold atoms as laboratory 
models of strongly-correlated electrons. In this respect 
the recent creation of fermion condensates (4] represents 
a significant step forward, as it allows the systematic ex- 
perimental investigation of one long-standing problem in 
that area, namely the crossover from BCS superfluidity 
of weakly-bound "Cooper pairs" to Bose-Einstein conden- 
sation of tightly-bound molecules jj| . 

One difficulty is that, although time-of-flight (TOF) 
experiments readily yield the momentum distribution of 
atoms in the condensate, n (fc), correlation functions are 
not simple to obtain. This is in contrast with electrons 
in solids, for which correlation functions can be obtained 
using neutron and X-ray scattering. In this context it 
has been proposed to use the pattern of atom shot noise 
in TOF images to deduce the density-density correlation 
function . Recent experiments have implemented these 
ideas to probe a Bose gas in the Mott insulating state 
and pairing correlations in a Fermi gas 0. In the case 
of a fermion superfluid, the density-density correlation 
function can be used to infer the superfluid wave function 
0. Here we point out that the latter quantity has very 
different properties in the BCS and Bose limits and show 
how to use it to identify the crossover region. 

We consider a Fermi gas where, as in a recent experi- 
ment , each atom can have two flavours of the "spin". 
For simplicity we assume translational invariance and 
work at zero temperature. The quantity of interest is 
the spin-flip density-density correlation function between 
two atoms with different momenta fik, tik': 



Figure 1: Qualitative dependence of the chemical potential, 
[i, and amplitude of the pairing potential, A, of a fermionic 
supefluid as a function of the strength of fermion-fermion at- 
traction, g. 



by averaging over the other magnitude |k'| and the two 
angles k, k'. This yields the |k|-dependence of the super- 
fluid wave function, <f>(k). Since (j> (k) = (c^^^^j it 

yields essential information on the pairing correlations in 
the ground state. It is an experiment of this latter type 
that we are considering here. 

Our key assumption is that the BCS ground state, used 
to derive JQ, adequately describes, at zero temperature, 
the correlations present in both the BCS and Bose lim- 
its 0. Using that ansatz, it is easy to show that the 
superfluid wave function is given by 



</>(k) = 



Al 



h 2 k 2 
2m 



A k 



(2) 



T (k; k') = (n T (k) n x (k')> - (n T (k)> (n x (k')) . 

Following Altman et al. [5] we evaluate this expression in 
the BCS ground state, to find 



r(k;k') = 0(k)<5 k ,_ k , 



(1) 



As pointed out in that reference the salient feature is the 
strong correlation between particles with opposite mo- 
menta, which can be observed in an experiment of the 
type described in 0, where the magnitudes |k|, |k'| and 
the angle k are integrated over and pairing correlations 
show up as a peak at k / = k + w. Additionally, it was 
suggested to study the dependence of T (k, — k) on |k| 



where A k is the "pairing potential" of BCS theory, m 
the mass of an atom and [i the chemical potential, /i 
and Ak have to be determined by solving the BCS self- 
consistency equations 0. Solutions can be found in 
the literature for a variety of models with s- and d- 
wave pairing, corresponding to specific choices of a pair- 
ing interaction potential V (k, k') — see, for example, 
Refs. Our emphasis here is on features that 

do not depend on the model under consideration. 

The qualitative behaviour of the pairing and chemical 
potentials, as functions of the strength g of the attractive 
interaction, is sketeched in Fig. For weak attraction, 
the chemical potential is independent of g and equal to 
the Fermi energy, E-p, and the amplitude of the pairing 
potential, A (where any k-dependence that does not vary 
with g has been factored out), grows like A ~ exp (— l/g). 
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As the attraction grows, A continues to grow, albeit more 
slowly, and the chemical potential eventually starts to 
decrease. This crossover regime sets in smoothly when 
A ~ fj,. In three dimensions, this happens near the value 
of g at which the interaction potential has a bound state 
and the scattering length diverges, g c . This lowering of 
the chemical potential corresponds to a shrinking Fermi 
surface. At some value g = go, the chemical potential 
goes below the bottom of the band: [i < 0. For g > go, 
the lowest-lying single-particle excitations occur at k = 
and there is no Fermi surface. Thus go is of more physical 
significance, from the point of view of the crossover, than 
9c 0. 

We expect the superfluid wave function, Eq. to re- 
flect the different nature of the BCS and Bose sides of the 
crossover. For simplicity we neglect the k-dependence of 
the pairing potential: Ak = A. Fig. |2] shows 4> (k) for 
different values of the ratio [i/A. The BCS limit cor- 
responds to /i > A. In this limit <f> (k) is maximum 
on the Fermi surface. As \x decreases, this surface be- 
comes smaller so the peak moves to lower values of k. 
For negative fx, (j> (k) becomes monotonic, with its max- 
imum at k = 0. This makes <j> (k) qualitatively different 
in the BCS regime (/i > 0) and the Bose-Einstein regime 
(M < 0). 

Because the maximum of cp (k) is at the locus of lowest- 
lying single-particle excitations, tracking the position of 
this maximum in k-space allows to monitor the grad- 
ual collapse of the Fermi surface along the BCS to Bose 
crossover. This can be achieved by comparing the maxi- 
mum value of <f> (fc) with its value at k = 0. The ratio is 
given by 



.(0) 




for fi > 
for fi < 



(3) 



This curve is plotted in the inset of Fig. [2 Note that it 
depends only on /u/A, not on the specific relation between 
A and fi. In that sense Eq. 10 is model- independent [T^ |. 

To further quantify the crossover in terms of (j> (k) , let 
us calculate its behaviour near k = 0. Differentiating (01 
we obtain d(f>/dk\k=o = and 



d 2 <j> 
dk 2 



fc=0 



2m (A2 + ^2)3/2' 



(4) 



This is plotted in Fig. We can see that the curvature 
of <fi (fc) at fc = changes sign when fi goes below the 
bottom of the band. Thus it provides a second test for 
the point on the phase diagram where the Fermi surface 
disappears. In fact for /i > Eqs. © and J3J, taken 
together, completely determine fi and A from measruable 
quantities. We note that the pairing potential can also be 
obtained directly from the large-fc behaviour of Eq. (2J; 
<f>(k) k (mA/h 2 ) k- 2 . 

To summarise, we have discussed, on the basis of the 
BCS ground state, the superfluid wave function of a 
fermionic condensate, <f>(k), which can be measured di- 
rectly in shot-noise correlation experiments. This has 
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Figure 2: The superfluid wave function for different values of 
At/A, as indicated. Inset: dependence of the maximum value 
of 4> (k) on At/A. 
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Figure 3: Contour plot of the curvature of the superfluid wave 
function at k = 0. The curves correspond to constant values 
of k%d 2 (j}/dk 2 \k=o, as indicated, k-p = (37r 2 n) 1//3 is the Fermi 
vector, given in terms of the average density of fermions, n, 
and £f = h 2 k 2 F /2m the Fermi energy. 



quite different behaviour in the BCS and Bose limits, 
changing qualitatively when the chemical potential goes 
below the bottom of the band. This point in the phase 
diagram marks the demise of the Fermi surface. Our 
main results are Eqs. 13141 . which describe the behaviour 
of the superfluid wave function along the BCS to Bose 
crossover. In particular, they allow to track the pair- 
ing and chemical potentials and to identify where on the 
phase diagram fi = 0. Our results are independent of 
any specific model of the interatomic interaction, relying 
only on the BCS ansatz for the ground state. 
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Note added. Since the author made this work avail- 
able in preprint form he has become aware of two un- 
published preprints discussing density correlations in 



the BCS to Bose crossover: W. Belzig, C. Schroll and 
C. Bruder, |cond-mat/0 412269; and B. Mihaila et al, 
cond-mat/0502110| 
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